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Abstract. Both the classical time-ordering and the Magnus expansion are well-known in the context 
of linear initial value problems. Motivated by the noncommutativity between time-ordering and time 
derivation, and related problems raised recently in statistical physics, we introduce a generalization of the 
£Nj . Magnus expansion. Whereas the classical expansion computes the logarithm of the evolution operator of 

a linear differential equation, our generalization addresses the same problem, including however directly a 
non-trivial initial condition. As a by-product we recover a variant of the time ordering operation, known 
as T*-ordering. Eventually, placing our results in the general context of Rota-Baxter algebras permits 
r*j ' us to present them in a more natural algebraic setting. It encompasses, for example, the case where one 

^ , considers linear difference equations instead of linear differential equations. 
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Introduction 

The time-ordered exponential Y(t) = Texp( f Q * A(s)ds^) solves the first order linear differential equation: 
X ' Y(t) = Y(t)A(t), 

with initial value Y(0) = 1. It is complemented by the solution written as a proper exponential of 
Magnus' expansion, Y(t) = exp (0(A)(t)). Therefore time-ordering and Magnus' expansion are naturally 
related. Indeed, the function 0(A) (t) is simply defined as the logarithm of the time-ordered exponential. 
It is characterized as the solution of a particular non-linear differential equation: 

and can be calculated recursively, which in turn is accompanied by formidable computational challenges 
due to intricate algebraic and combinatorial structures. 

Linear initial value problems are abundant in the mathematical sciences. And therefore the time- or 
T-ordering operation plays an essential role in many fields, for instance in quantum theory and statistical 
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physics. In concrete applications one often needs to calculate time derivations of time-ordered prod- 
ucts. However, the noncommutativity of the T-ordering operation and time derivation causes obstacles. 
Our work started out with the aim to understand a simple though surprising paradox caused by the 
aforementioned noncommutativity. It turns out that physicists introduced a variant of the time-ordering 
operation, denoted T*, which precisely encodes the lack of commutativity between the time-ordering 
prescription and time derivations. As part of this work we explore the relation between this modified 
time-ordering operation and the Magnus expansion (Theorem[T]). This lead to a generalization of Magnus' 
result, allowing to include the initial value into the Magnus expansion (Theorem [2]) . 

One may consider the time-ordered exponential as the formal solution of the linear integral equation: 



corresponding to the initial value problem. It recently became clear |12] that such integral equations, 
seen as linear fixpoint equations, are just a particular example of a more general class of abstract linear 
fixpoint equations: 



defined in associative Rota-Baxter algebras (Theorem[5]). We therefore embedded the generalized Magnus 
expansion into the context of such algebras, by showing in detail a theorem that specializes to the 
particular case (Theorem [7]) . 

Let us remark that both time-ordering and Magnus' expansion have triggered much attention and 
progress in applied mathematics and physics. See for instance [6j [171 [20J, (26J [29]. Recently, however, 
Magnus' work has also been explored successfully from a more algebraic-combinatorial perspective using 
operads, rooted trees, Hopf algebras, pre-Lie algebras, and the theory of noncommutative symmetric 
functions El E21 ESI EU ■ 

Finally, as a word of warning, we should mention that we deliberately ignore any convergence issues. 
This work addresses mainly formal algebraic and related combinatorial aspects. 

The paper is organized as follows. In the first section we introduce the objects we are going to consider 
in this work, i.e. the time-ordering map and the Magnus expansion. In section [2] we introduce a variant of 
the time-ordering operation, which reflects the non-commutativity of time-ordering and time derivation. 
The connection to the Magnus expansion is detailed. Section [3] recalls briefly the notion of Rota-Baxter 
algebras, and we then embed the foregoing results into this context. We also show a mild generalization of 
Atkinson's factorization theorem for Rota-Baxter algebras. The notion of -R-center is introduced, which 
plays an important role. Finally, in section [4] we present finite differences as a simple example. 



1.1. The T-ordering map. Linear initial value problems (I VP) play a major role in many parts of 
applied mathematics and theoretical physics. They consist of an ordinary linear differential equation for 
an unknown function Y(t): 



together with a specified value of this function at a given point in time, e.g. say for t = 0, Y(0) = Yq, which 
we call the initial value. The functions Y(t) and A(t) being for instance matrix valued and differentiable. 

Further simplifying the problem, e.g. by restricting to scalar valued functions, the solution of (p} is 
given straightforwardly in terms of the exponential map: 




Y = l + R{AY) 



1. Time-ordering and the Magnus expansion 



(1) 




(2) 
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Returning to the more challenging case, when noncommutativity prevails, we see that the result is far 
from simple. Recall that the formal solution of (TJJ is given by the series of iterated integrals: 

(3) Y(t) = Y (l+ [ A( Sl )d Sl + f f 1 A(s 2 )ds 2 A( Sl )d Sl + •■ 

V Jo Jo Jo 

where 1 would be the identity matrix in the context of matrix valued differentiable functions. The right 
hand side, known as Dyson-Chen or T(ime)-ordered series, follows forthrightly from the solution of the 
integral equation: 

(4) Y (t) = Y + [ t Y(s)A(s)ds 

Jo 

corresponding to ([I]). Note that in the commutative case series (|3|) simply follows from expanding the 
exponential map ((2|) and making use of the integration by parts rule: 

rt rt rt rs rt ru 

(5) / f(s)ds / g(u)du= / / f(s)g(u)duds + / f(s)g(u)dsdu. 
Jo Jo Jo Jo Jo Jo 

Let us recall the notion of time- or T-ordering. Time ordering is defined in terms of the T-map: 



U(si)V(s 2 ), si < s 2 
V(s 2 )U(si), s 2 < si. 



J[U( Sl )V(s 2 )] 

Written more concisely it involves Heaviside step functions |15j : 

T[U( Sl )V{s 2 )] := 0(s 2 - Sl )U{ Sl )V(s 2 ) + 8(*i - s 2 )V(s 2 )U( Sl ), 
with a generalization to higher products of operators: 

T[C7i(ai) ■■■£/„(««)] := ®nUa(i)(sa(i))---U ain) (s a{n) ), 

where 0^ := IliLi ®( s o-(i) — s ct(2)) " " " ®( s o-(n-l) — s a(ri))- Observe that for any permutation a € S n , 
T[J7i(si) • • • U n (s n )] = TfC/^n^s^i)) • • • U a (n)(. s iT(n))]- Following the convention that the T-map operates 
before integration, the solution of the IVP ([TJ writes: 

Y(t) = Y Texp ( f A(s)ds) := Y V 1 [ ■■■ [ J[A( 8l ) ■ ■ ■ A(s n )]d Sl ■ ■ ■ ds n . 
\Jo J n > n[ Jo Jo 

Each n-fold T-product inside the integral produces n! strictly iterated integrals, and therefore the series 
coincides with ©. 

As a remark we would like to mention that our work started out with the aim to understand a 
simple though surprising paradox caused by the noncommutativity of the T-ordering operation and time 
derivation. See [3] for an interesting example where this plays a decisive role, and |27j for a case where this 
paradox caused some confusion. We assume nonnegative t, and At := A(t) may be a finite-dimensional 
matrix valued function, or more generally a linear operator valued map. The following identity: 



exp (A t - A Q ) = exp ^ 



au- 



dit 



may be considered as evident. However, its T-ordering cousin, which would result from a naive action of 
the T-ordering operation on both sides of the foregoing equation: 

/ ft dA 

(6) Texp (A t — Aq) = Texp I / du 



du 
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does not hold. Indeed, let us look at the second order term of each side separately: 



T[(A t - A)) 2 ] = A A - 2A A t + A t A t , 



whereas: 



T 



dA 



u n2 



du 
t 



dv- 



+ 







du 







du dv 

^ v ^"^ u dA v ^ 

du dv 



du / dv 



dv du 
du^(A t -A u ) 



o 







du 



Hence, we have verified that: 



= 2{A t - Ao)A t - 2 J dJ^A u . 

J[(A t -A f)^j\([ 
. Jo 



, dA u . 2 



In fact, it is clear that © does not hold due to the lack of commutativity of the T-ordering with time 
derivation, which, on the other hand, is caused by the Heaviside step functions involved in the definition 
of the T-ordering map. However, observe that the difference reads: 



T[(A t -A y 



o 



'du^f 



du 



-A t A t + A A + 2 



du 



o 



dA u 
du 



A u 



-A t A t + A A + 



du 



u A u + A u ^ | + 



du 



du 



o 



* du ( d 4^A u -A u ^ 
du 



-A t A t + A A + A t A t - A A + 



du 



dAu 
du 



A u A u 



du 
dAu 
du 





'dA u ' 


/ du 




Jo 


du 



In the sequel we will further remark on this observation, and explore in more detail the object, denoted 
^K^f)' sucn that the following equality holds: 

/ f 4 

(7) 



Texp (A t - Aq) = Texp (J du *(^)) • 



1.2. The Magnus expansion. It is clear that the proper exponential solution of (HJ changes drastically 
in the light of the noncommutative character of the problem. Wilhelm Magnus proposed in his seminal 
1954 paper [19] a particular differential equation for the logarithm of the series ([3]) of iterated integrals 
(with initial value Yq = 1), which we denote by the function £l(A)(t): 



B 



n(A)(t) = A(t) + Y,(-lT^aSp 



n>0 



n \ JJn(A)(a)efa 



(A(s)) 



-ad 



n(A) 



g-a d n(A) 



with tt(A)(0) = 0. The B n are the Bernoulli numbers -Bo — 1) B± ■ 
for k > 1, such that the solution of the IVP ([T|) (with Yq = 1) is given by: 



2' 



Bo 



and B 2 k+i 



(8) 



Y(t) = exp / n(A)(s)ds 



Here, as usual, the n-fold iterated Lie bracket is denoted by ad^a\b) := [a, [a, ■ ■ ■ [a, &]]■■■]. Recall that: 



exp(— z) 



and 



exp(— z) — 1 



n>0 



exp(—sz)ds. 
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The reader is encouraged to consult [6] for an authoritative review on the Magnus expansion and its 
ramifications. References |20[ [2"6[ [29] are classical sources regarding Magnus' work. Let us write down 
the first few terms of what is called Magnus' series, f2(AA)(s) = Xm>i ^ n ^n(A)(s), following from Picard 
iteration to solve the above recursion: 



j t n(\A)(t) 



XA(t) + A 2 - 



1 



A(s)ds,A(t) 



A(u)du,A(s) ds,A(t) 



+ A 3 — 
12 



A(s)ds, A(u)du,A(t) 



+ 



where we introduced a parameter A for later convenience. In the light of the IVP (JTJ) and its two 
solutions, i.e. the proper exponential solution using Magnus' expansion, and the T-ordered exponential, 
we find the identity (for Yq = 1): 



(9) 



Texp 



A(s)ds 



exp 



n(A)(s)ds . 



The terms beyond order one in £l(\A)(s) = A J * A(s)ds+^2 n>2 \ n Q n (A)(s) consist of iterated Lie brackets 

and integrals. For each n > 1 they solely serve to rewrite higher powers ( J" * A(s)ds) n as strictly iterated 
integral of order re, modulo a factor re!. For instance, verify that: 



A(s)ds + 



1 



A(u)du,A(s) 



ds 



o Jo 



A(u)du A(s)ds 



We remark that recently, Magnus' work has been successfully explored from a more algebraic-combinatorial 
perspective using the theory of operads, rooted trees, pre-Lie algebras, Rota-Baxter algebras, and the 
theory of noncommutative symmetric functions [HI [TTj fT2[ \16\ \TT\ [21] . 



2. T-ORDERING, T*-ORDERING AND TIME DERIVATIONS 

2.1. T-ordering and time derivations. Let us return to the erroneous identity 
simplify the presentation, we sometimes write At for A{t). 

Theorem 1. The solution to the equation: 

Texp (A t - Aq) = Texp du *(^)) 

is explicitly given by: 



Recall that to 



(10) 



n=l m=0 



du 



{Ay) 



dAu 
du 



{An) 



n—l—m 



-od Au _ 1 / dAu 



-ad An 



du 



Observe that if A u commutes with its derivation, [^^S^-u] = 0> then, as expected equation (|10p 
reduces to $(^) = ^. 

In fact, equation (|10[) follows from exploring briefly what went wrong in the first place with identity 

([6]). Indeed, recall that Xt = Texp ( duB^j for a family of operator valued maps Bt is the solution to 
the simple IVP ±X t = X t B t , X = 1. Hence: 



d_ 

di 



Texp 



' dA u 
du— — 
du 



Texp 



du 



dA u \ dA 



du j dt 



On the other hand it is clear from its definition, that: 

Texp (A t - Aq) = exp (-A ) exp (A t ) . 
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Therefore: 

(11) ^-Texp (A t - A ) = exp (-A ) exp (A t ) = exp (-A ) exp (A t ) exp (-A t ) ^- exp (A t ) 
at at at 

= (Texp(^ i -A ))exp(-^ t )— exp (At) , 

which again shows that ([6]) can not be correct. 
However, the last equality may be used to derive ([TO]) . Indeed (fTT|) implies: 

Texp (A t - A ) = Texp ( / duexp (-A u ) exp (A u 

\Jo au 

Using Duhamel's formula: 

d f 1 dA 

(12) exp (-Ay) — exp (A u ) = / dv exp (-vA u ) — ^ exp (uA u ) , 

J 

and recalling that: 

-l)'" 1 .,, e~ ad °- - 1 



/ di/exp(-i/a)6exp(i/a) = V ^ ad**" 1 ^) = ; — (6), 



yields: 



e -ad Au _ 1 , dAu 



-adA u \ du 



(13) exp (—Ao) exp (At) = Texp (A t — Aq) = Texp (^j du 

This implies (TO]) because ad^ (b) = £™~ n (-l) n C^a n ba m - n . Conversely, if we set B t := $(^f) and 
:= exp (Aj), we arrive at the following first generalization of the classical result of Magnus ©. 

Theorem 2. The IVP: X t = X t B t , X = exp(a) is solved by: 

X t = exp (® + J & a (B)(s)ds^J , 

with: 

Cl a (B)(u) := -j ^ B u . 



2.2. T*-ordering and the Magnus formula. Recall that the modified time- or T*-ordering is defined 
precisely so as to encode appropriately the lack of commutativity between the T-ordering and time 
derivations. In [5] the interested reader can find motivations for the introduction of this particular 
variant of the time ordering operator coming from statistical physics, as well as further bibliographical 
references on the subject. 
For I(t) := d( ^ d ^ , we define: 

T* [1(h) ■ ■ ■ i(t n )] ■.= A • • ■ -^-JiQih) ■ ■ ■ Q(t n )]. 

For any permutation a £ S n : 

T*[J(t!) •■■/(*„)] =T*[/(i ff(1) ).../(i ff(n) )], 

or, equivalently, we have an invariance under the right action of the symmetric group, S n , in n elements, 
T* = T* o a. In particular, because of this property, T* is entirely characterized by its action on 
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time-ordered products I(t\) ■ ■ ■ I(t n ), t\ < ... < t n , or, equivalently on products I (ti) mi ■ ■ ■ I (t n ) r 
h < ■ ■ ■ < t n . 

Therefore we get, by setting I(t) = A(t) and Q(t) = A(t) - A(0): 

t t 

T*ex V yy u A u "j = ^^J...jfldt^--~J[(A(t 1 )-A(0))...(A(t n )-A(0))) 

n Q i=l 

n 

= E E ^ (f) (-V k A(0) k A(tr- k = exp(-A(0)) «p(A(t)). 

n k<n ' ^ ' 

Equation (fl~3j) yields: 

(14) T*exp (J duA^j = Texp (J du & 1 (A u ) 

Let us look at the degree two terms in A and A on both sides of this equality. We find: 

— / / T*[A Ul A U2 )duidu2 = / / T[A Ul A U2 ]duidu 2 + — / ad Aui {A Ul )dui. 

' J ** * « " * ^ 

Simply rewriting the last term on the right hand side: 

— / adA (A Ul )dui = — / ad A {A U2 )5(ui - u 2 )duidu 2 
z - Jo l - Jo Jo 

we get finally, for u\ different from u 2 : 

T \A U1 A U2 \ = T[A U1 A U2 \. 

On the diagonal, i.e. for u\ = u 2 we obtain: 

T*[A Ul A U2 ] = ad Aui (A U2 )5(ui -u 2 ). 

In general, by identifying homogeneous components on each side of the foregoing equation, and taking 
into account the permutational invariance of T* (so that e.g. T*LA Ul ^4 U2 yl Ul ] = T*[A^ A Wa ]), we finally 
arrive at: 

r, \ . ad^ . ad ( ; fc " 1} . 1 

#1 • • ■ 8k, 



n\\nx,...,n k J 



ad^ ad { : k ~ l) 

(-ir-zn— ■ ■ ■ (-ir— ^(A Ck ) 



for arbitrary < c\ < ■ ■ ■ < c k < t, and n±, . . . , n k , such that n = n\ + ■ ■ ■ + n k , u\ = ■ ■ ■ = u m = 

ci, . . . ,«niH hnfc_i+i = • • • = u ni -\ \ rnk _ 1+nk = Cfc. Recall that the multinomial coefficient on the left 

hand side stands for n\{n\\ • • • n k \)~ l ■ Here we wrote Si for the product of delta functions 5(u ni ^ i-nj_i+i — 

u ni+ ... +n ._ 1+2 ) ■ ■ ■ 5(u ni+ ... +rii _ 1+ni _ 1 -u ni+ ... +rH _ 1+ni ) encoding the change from a n-dimensional integral 
to a fc-dimensional integral from the left to the right hand side of equation (|14p . Note that we have 
recovered [H eq. 2.11] (up to signs, due to different conventions regarding time-orderings) : 

Corollary 3. We have, with our previous notation: 



(15) T*[A U1 ■ ■ ■ A Un ] = (-irJ[ad^'(A Cl ) ■ • • a^^A^h -6, 
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3. Linear recursions in Rota-Baxter algebras 

We take now a broader perspective on linear IVP by replacing ^ in terms of general linear operator 
fixpoint equations. A natural setting to do so is provided by associative Rota-Baxter fc-algebra, which we 
briefly recall in this section. For more details we refer the reader to [3"1 [5"! 110 1 114 1 124] and references therein. 
See also [H [2] for interesting aspects. We assume the underlying base field k to be of characteristic zero, 
e.g. k = E or C. 

Motivated by Frank Spitzer's seminal work [25] in probability theory, the mathematician G. Baxter [5], 
soon after followed by J. F. C. Kingman, W. Vogel, F. Atkinson, P. Cartier, Gian-Carlo Rota and 
others [31 [71 [181 ESI [ZS [23l [23] , explored a more general approach to the above linear IVP (p} and the 
corresponding linear integral equation ^) - in the case of lo = 1 ~i by suggesting to work with the linear 
fixpoint equation: 

(16) Y = 1 + XR(Ya), 

in a commutative unital /c-algebra A, a € A, respectively the filtered and complete algebra *4[[A]]. The 
linear map R on A is supposed to satisfy the Rota-Baxter relation of scalar weight 8 £ k: 

(17) R(x)R(y) = R{R{x)y + xR{y)) + 9R{xy). 

We call A a commutative Rota-Baxter k- algebra of weight £ k with unit 1. Observe that the map 
R := —6 id — R is also a Rota-Baxter map of weight 8. Moreover, from identity (|17p it follows that the 
images of R and R are subalgebras in A. 

Baxter gave a solution for (|16p by proving Spitzer's classical identity in a commutative unital Rota- 
Baxter /c-algebra of weight 9 £ k, that is, he showed that the formal solution of (|16p is given by: 



(18) Y = 1 + R(R(- ■ ■ R(R (a)a)a ■ ■ ■ )a) =exp(-R ^° S ^ g °M 



n=l 

ra-times 



Remark 4. Note that besides (j!6p . Baxter, Atkinson and Vogel considered more general linear fixpoint 
equations and their solutions in a commutative Rota-Baxter /c-algebra A. In [31 [5] the equations E = 
a + \R(bE), F = R(a) + XR(bF), with analog equations for R, and G = a + XR(bG) + XR(cG) are solved 
for a, b, c £ A. In [28j Vogel considered the equation H = a + XbR(H) + cXR(H), and gave a solution for 
a,b,c £ *4. The generalization of these equations and their solutions to non-commutative Rota-Baxter 
algebra was given in |121 [13] . 

Let us return to (|17p . which one may think of as a generalized integration by parts identity. Indeed, the 
Riemann integral, which satisfies the usual integration by parts rule ([5]), corresponds to a Rota-Baxter 
operator of weight zero. Observe that in the limit 9 — > we obtain: 



1 lo g (l-0a) = -J2d n - 1 - " 

n>0 



The right hand side of the last equality in ([18]) reduces to the simple exponential solution Y = exp(R(a)). 
In the context of the Riemann integral, this corresponds to the classical case ([2]) of the IVP ([1]), with 
Y = l. 

The extra term on the right hand side of (j!7p becomes necessary, for instance, when we replace the 
Riemann integral by a Riemann-type summation operator: 

(19) £(/)(*) + 

n>0 
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on a suitable class of functions. Indeed, one verifies that the map S satisfies the weight 9 = 1 Rota-Baxter 
relation: 

S(f)S(g) = S(S(f) g)+S{f S(g)) + S(fg) . 
More generally for finite Riemann sums: 

[x/9]-l 

(20) R 9 (f)(x):= Y e ^ n0 )- 

n=0 

we find that Rg satisfies the weight 6 Rota-Baxter relation: 

Re(f)Re(g)(x) = Rg{Rg{f)g){x) + Rg{fRg{g)) (x) + 9R e (fg)(x). 

Beside fluctuation theory in probability the Rota-Baxter relation recently played a crucial role in Connes- 
Kreimer's Hopf algebraic approach to perturbative renormalization [9]. Here, Rota-Baxter algebras enter 
through projectors, which satisfy (flTl) for 9 = —1. A paradigm is provided by the algebra A = Cfe" 1 , e]] = 
e _1 C[e _1 ] © C[[e]] of Laurent series X^L-fc a n e ™ (with finite pole part). The projector, called minimal 
subtraction scheme map, used in renormalization is defined by keeping the pole part: 

oo —1 oo —1 

Rms( Y a ' n£n ) Y a n £ n + Y a n^ n ) ■= Y an£n - 

n=—k n=—k n=0 n=—k 

It is a Rota-Baxter map of weight 9 = —1. In general, assume the £>algebra A decomposes, A = Ai®A2, 
and let R : A — > A be defined by R(a±, 02) = ai, then R 2 = R. Hence: 

R{a)b + aR{b) - ab = R(a)R(b) - {id - R){a){id - R){b) 

such that applying R yields i?(i?(a)6 + aR(b) — ab) = R(a)R(b). We refer the reader to Q3] for more 
details including examples and applications. 

In |12t [13] the Spitzer identity has been generalized to non-commutative Rota-Baxter algebras: 

Theorem 5. |12j Let (A, R) be a unital Rota-Baxter k-algebra of weight 9 € k. Let Q' := Q'(Xa), a € A, 
be the element o/A^4[[A]] such that: 

Y = exp( J R(ft / )), 

where Y = Y(a) is the solution of the linear fixpoint equation Y = 1 + XR(Ya). This element obeys the 
following recursive equation: 

(21) n'(Xa) = Z^m±I^ {Xa ) = y (_i r ^ al^\\a) 

e 1 m>0 

with B[ the Bernoulli numbers, rg a (b) := 9ba, and ad a {b) := adn( a \(b) — rg a (b). 

Remark 6. Observe that the right multiplication rg^i map accounts for the extra term in (|17p . Let 
us emphasize that in [12] we formulated this theorem using the generic pre-Lie algebra structure under- 
lying any associative Rota-Baxter algebra. Indeed, the above generalized adjoint operation, ad a {b) = 
a dR( a )(b) — rg a (b) = [R(a),b] — 9ba, defines a pre-Lie product on A. However, we refrain from using the 
pre-Lie picture for reasons to become clear in the sequel. 

We introduce now the notion of the i?-center in a Rota-Baxter algebra A of weight 9, defined as the 
following subset of elements in A: 

Zr(A) := {x € A I xR(a) = R(xa) and R{a)x = R(ax), Va € ^4}. 

Observe that the i?-center forms a subalgebra of A. Indeed, for n, v € Zr(A), we have uvR(a) = 
uR{va) = R(uva), and [R(a),x] = R([a,x]). Notice the identity R(x) = xR(l), for any x € Zr(A). In 
the context of matrix valued differentiable function with the Riemann integral as weight zero Rota-Baxter 
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map, the fZ-center includes amongst others the constant matrices. Another useful example are upper 
(or lower) triangular matrices, say, with unit diagonal, and entries in an arbitrary Rota-Baxter algebra 
A of weight 9. They form a noncommutative Rota-Baxter algebra -M4 of weight 9 with the natural 
Rota-Baxter map. In this case the triangular matrices with entries strictly in k are in the i?-center. 

Now we would like to prove that the combinatorial structure underlying the T*-operation, that is the 
solution of the IVP provided by Theorem [21 extends to Rota-Baxter algebras. 

Theorem 7. Let (A,R) be a unital Rota-Baxter k-algebra of weight 9 £ k. Let a be an element in its 
R-center, and Q' a := fL' a (Xa), a £ A, be the element of A[[X\] such that: 

X = exp(a + R(Q' a )), 

where X is the solution of the linear fixpoint equation: 

X = exp(a) + XR(Xa). 

The element Q,' Q obeys the following recursive equation: 

(22) W a (Xa) = (Xa) = V(-l)"^ ad^, (Xa), 

m>0 

where ad^Jb) := ad a+R[n , a) - r 6 n' a - 

Observe that this theorem is a generalization of Theorems [2] and [5j Returning to Remark [61 we see 
here that this generalization involves the operation ad a ^ := ad a+R ^- ) —rgb, which does not define a pre-Lie 
algebra product on the Rota-Baxter fc-algebra A. 

Proof. We follow partly [12]. First, note that we can decompose ad a+R ^i^ = £ a +R(n' a ) ~ r a+R(n' a )i where 
ta(b) '■= oh. It is clear that both operations commute, that is, i a fb( c ) = aco = r b^a(c). Now, let: 

Xa = e —— -(n' a ) = / dse- siad <*+w a )- r °<\n' a ) 

1 pi 

d se ~ s( ~ e c l +R(si' a )- r a +R(n> a )- r e(i' a )(£}' a ) = / dse~ sea + R< - n '^e sra + R ^ +sresi 'a (£l' a ). 
Jo 

Now, we multiply by X = exp(a + R(fl' a )) from the left: 

Xa = dse il ~ s)ia + R (<'>e s{ra + R< -<'> +re < ) (n' a ) 
Jo 

= f ds E {± ^(^R(n>J P (r« + R(nij + W 

P,<?>0 

p,q>o 

= E^y E (<x + R(n' a )y(n> a )(a-R(n> a )y, 



n>0 p+q=n—l 



where we used that —R = 9id+R. Hence, we would like to show that X — exp(a) = XR(Xa). Expanding: 

X - exp(a) = ^ 1 (( a + R(W a )) n - a n ) , 



n! 

n>0 
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lefts us with the goal to prove order by order the general identity: 

n>0 n>0 ' p+q=n— 1 

where a € Zr(A) and /3 is an arbitary element in A. 

The case n = 1 is obviously true. Let us look at n = 2. Then we have on the left hand side: 

{a + R(P)) 2 -a 2 = aR{P) + R{p)a + R(P)R(P) 
= R{a(3) + R(p a ) + R(p)R(p). 

In the last line we used that a is in the i?-center Zr{A). On the right hand side we find: 

R{{a + R{p))p + P{a-R{p))) = R(a(3 + R(p)p + /3a + pR(J3) + 0/3/3) 

= R(aP + Pa)+R(P)R(P). 

In the last line we used the Rota-Baxter relation (|17p . We show now the case n + 1, assuming the 
identity holds up to order n. This yields: 

(a + R(P)) n+1 - a n+l = (a + R(P)) n {a + i?(/3)) - a n a 
= (a + R(P)) n R(P) + ((a + R(P)) n -a n )a 
= £ i?((a + i?(/3)) p (/3)(a-i?(/3)) g ),R(/3)+a ri ii(/3) 

p+q=n~ 1 

+ £ ^((a + ^(/3)) p (/3)(a-fi(/?)) 9 a) 

p+q=n— 1 

= £ i?((« + ^)f(/3)(a-E(/3))^(/3)) 

p+q=n~ 1 

+ ^ R(R{(a + R(P)Y(P)(<*-R(P)) q P)) 

p+q=n—l 

+ Y R{(a + R(P)) p (P)(a-R(p)yep) 

p+q=n—l 

+R{a n p)+ Y R{(a + R(P)) p (P)(a- R{P)) q a) 

p+q=n—l 

= y R{(a+Rm p ma-R(p)) q+1 ) 

p+q=n— 1 

+ £ i?(i?((a + ii(/3))P(/3)(a-i?(/3)) 9 )/3) + fl(a n /3) 

p+q=n— 1 

= £ i?((« + i?(^(/3)(a-E(/3))«+ 1 ) 

p+g=n— 1 

+R((a + R{P)) n P - a n p) + R{a n p) 

= Y R(^ + RW)) p W)(a-R(p)) q+l )+R(^ + R(P)) n P) 

p+q=n~l 

= Y R{^ + R{P)) p m^-R(P)) q ) 

p+q=n 

where we recall that a — R{P) = a + OP + R(P). □ 
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Note that in the context of the weight zero RB algebra defined in terms of the Riemann integral, 
the generalization following from the theorem allows to absorb the initial value in ([1]) directly into the 
Magnus' expansion. 

Remark 8. Recall Remark 0] and the fact that any Rota-Baxter algebra has two Rota-Baxter operators, 
R and R. In [U [5] Baxter and Atkinson showed that the solution to the recursions: 

E = a + \R{bE) and F = R(a) + XR(bF) 

are given by F = R(aX)Y and E = a + R(abX)Y , respectively, where X, Y are solutions to: 

X = 1 + XR(Xb) resp. Y = 1 + XR(bY). 

Together with the solutions X and Y for the last two fixpoint equations comes the so-called Atkinson 
factorization, which generalizes straightforwardly in the context of the above theorem. 

Proposition 9. Let (A, R) be a unital Rota-Baxter k-algebra of weight € k. Let a, f3 be in its R-center 
Zr{A). Then we have a factorization: 

1-Xa = X~ 1 exp( 7 )y~ 1 , 

where: 

X = exp(a) + XR(Xa) Y = exp(/5) + XR(aY), 

and exp(7) = exp(a) exp(/3), 7 G Zr(A). 

Proof. The proof reduces to a simple verification. 

XY = exp( 7 ) + XR(Xa) exp(/3) + A ex.p(a)R{aY) + X 2 R{Xa)R{aY) 
= exp( 7 ) + XR(Xa exp(/3)) + XR(exp(a)aY) + X 2 R(Xa)R(aY) 
= exp( 7 ) + A«(Xa(exp(/3) + XR(aY))) + A,R((exp(a) + XR{Xa))aY) 
= exp( 7 ) + XXaY, 

from which the result follows. □ 

4. Example: Finite differences 

As an example of an interesting fixpoint equation involving a non-zero weight Rota-Baxter structure, 
we will now approach the particular case of a linear finite difference initial value problem. For this, we 
define the finite time-difference operator A such that for all functions / defined on N: 

A(/)(n) := /(n+l)-/(n). 

Observe that this map satisfies a generalized Leibniz rule: 

A(/s)(n) = (/j)(n + l)-(/s)(n) = /(fl + l)s(fi + l)-/(n)j(n) 
= f(n)A(g)(n) + A(f)(n)g(n) + A(f)(n)A(g)(n) 

n-l 

Recall the finite Riemann sum map (|20|) . R{X) n := X^. The latter satisfies the Rota-Baxter relation 

k=0 

of weight 6. We put the weight equal to one, R := R\. Then: 

m— 1 m— 1 

fl(A(/))(m) = j>(/)(n) = ^(/(n + l)-/(n)) 

n=0 n=l 

m m—l 

n=l n=0 
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= /(m)-/(0). 
We are interested in the finite difference initial value problem: 

/nn\ AA n X n B n 

1 ' X = e W (A ) 

where X n := X(n) and B n ^ — 1. Applying i? leads to a generalized Atkinson equation: 

X = exp(A)) +R(XB), 

which yields: 

n-l 

X n = exp{A ) + J2 x kB k - 

k=0 

Observe that the solution can be written as a finite product, as well as, using Theorem^ as an exponential: 

n-l 

X n = exp (A)) J] ( J + B k) = ex P W(»))) • 

fc=0 

Now, by applying Theorem with a = A(0) and fi^ = AA (recall that X = exp(^4) and R(AA) = 
A - A(0)), we get: 

Theorem 10. With our previous notation, for the finite difference IVP 123\) we get: 

y —ZA h + r A k+1 J 

where as before £x is the left multiplication by X and rx is the right multiplication by X. 
In particular: 

(24) exp (-A ) exp (A N ) = ][| (l + ^ " " ) t A ^) 

Note that these formulas can also be obtained from a straightforward calculation. We include a direct 
proof for the sake of transparency. The starting point is a discrete analog of Duhamel's formula (|12p : 

Lemma 11. We have: 

(25) exp (-A k ) A [exp (A k )} = f dv exp (-vA k ) (AA k ) exp (vA k+1 ) . 

Jo 

The proof follows by taking the integral from to 1 of the following equality: 

^- [exp (-vA k ) A [exp (yA k )]] = - exp (-vA k ) A k A [exp (vA k )] + exp (-vA k ) A [A k exp (iM fe )] 
dv 

= - exp (-zM fc ) A k A [exp (zM fc )] + exp (-vA k ) A k A [exp (iM fc )] 

+ exp (-vA k ) A [A k ] exp (vA k ) + exp (-i/A fc ) A [A k ] A [exp (iM fc )] 
= exp (—vA k ) A [A k ] exp (uA k+1 ) . 

Now we verify identity (|24p . The first thing to remark is that: 

JV-l 

exp (-A ) exp (,4jv) = JJ (1 + exp (-A fc ) A [exp (A k )}) 

k=0 
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follows when we note that 1 + exp (— Aj.) A [exp (A/.)] = exp (— A^j exp (Ak+i) ■ Then, we first apply 
Duhamel's formula ([25]) . and second expand exp(— vA^) and exp (vAk+i) , and eventually we do the 
explicit integration of v. This yields: 



N-l / oo „ N-l 



(26) exp (-Ao) exp (An) = II 1 + E E W A M (^W*" 1 " 



iV! 

fc=0 \ AT=1 n=0 



Finally, we note that: 

N-l 
n=0 



such that: 



oo iV-l oo /-i-vJV-l 

E a77 E (- l T C N-i (A k ) n A[A k ] (Ak+xf- 1 - 71 = EM^T-^-^J^A^] 



jV=l ' n=0 N=l 



exp(-/ Afc +r Ak+l ) - Id 
-U k + r Ak+1 



A [A k ] . 



This relation, together with (I26p concludes the proof of ([24 
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